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Summary 


This  report  presents  an  application  of  tedaniques  developed 
under  this  contract  to  use  regional  seismic  network  data  to  define  the  varia¬ 
tion  of  Lg  amplitude  with  distance  and  frequency. 

The  distinguishing  feature  of  this  work  is  that  we  do  not  use  a  simple 
relation  for  Lg-  or  S-wave  amplitude  variation.  Rather  we  use  a  very  general 
form 

a(/,  r) = , 

where  the  geometrical  term  g{r)  is  not  a  simple  functional  relation  such  as 
r  -  71.  To  permit  the  observed  data  themselves  to  define  the  empirical  form, 
typically  1000  -  2000  wave  forms  are  used  in  the  distance  range  of  interest. 
Regression  is  used  to  define  a  distance  dependence  which  is  then  interpreted 
in  a  second  stage  in  terms  of  geometrical  spreading  and  Q.  It  is  at  the  second 
stage  that  a  prediction  model  is  developed  using  random  vibration  theory  as  a 
tool  that  connects  Fourier  spectra  and  signal  duration  to  peak  time-domain 
amplitudes. 

The  surprising  feature  that  we  have  noticed  in  analyzing  data  from  the 
southern  Great  Basin  surroimding  NTS,  southern  Cahfomia,  the  Pacific 
Northwest,  the  central  United  States,  Germany  and  Italy,  is  the  difference  in 
the  g{r)  at  distances  greater  than  40  km  fi:om  the  source  in  addition  to  the 
expected  differences  of  Qif).  These  variations  must  be  the  result  of  regional 
differences  in  crustal  structure.  Our  data  depend  on  the  structure  but  are 
insufficient  to  define  the  structure. 

The  techniques  used  here  are  applicable  to  those  parts  of  the  world  with 
regional  seismic  network  observations  of  local  events.  While  this  constraint 
precludes  the  universal  use  of  our  techniques,  they  are  useful  in  some  regions 
of  interest. 
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1.  Introduction 


Ground  motion  attenuation  with  distance  and  the  variation  of  excitation 
with  magnitude  are  parameterized  using  three-component,  0.25-5.0  Hz  earth¬ 
quake  ground  motions  in  the  distance  range  of  15  -  500  km  for  southern  Cali- 
forma.  The  data  set  consists  of  820  three-component  TERRAscope  recordings 
from  140  earthquakes,  recorded  at  17  stations,  with  moment  magnitudes 
between  3.1  and  6.7.  Regression  analysis  uses  a  simple  model  to  relate  the 
logarithm  of  measured  ground  motion  to  excitation,  site,  and  propagation 
effects.  Regression  is  performed  on  Fourier  velocity  spectra  and  peak  veloci¬ 
ties  in  selected  narrow  band  pass  filtered  frequency  ranges.  Regression 
results  for  Fourier  ampHtude  spectra  and  peak  velocities  are  used  to  define  a 
piecewise  continuous  geometrical  spreading  function,  frequency  dependent 
Qif))  and  a  distance  dependent  duration  that  can  be  used  with  random  vibra¬ 
tion  theory  (RVT)  or  stochastic  simulations  to  predict  other  characteristics  of 
the  ground  motion. 

The  duration  results  indicate  that  station  terms  are  required  and  that  both 
the  degree  of  distance  dependence  and  scatter  decrease  with  increasing  fre¬ 
quency.  The  ratio  of  horizontal  to  vertical  component  site  terms  is  about  V2 
for  all  frequencies.  However,  this  ratio  is  nearer  to  unity  for  rock  sites  and 
higher  for  soil  sites. 

Modeling  indicates  that  the  Fourier  velocity  spectra  are  best  fit  by  bilinear 
geometrical  spreading  exponents  of  -1  for  r  <  40  km  and  -1/2  for  r  >  40  Irm 
The  frequency  dependent  quality  factor  is  Q  (f)  =  180  for  each  of  the 
three  components  and  also  for  the  combined  three-component  data  sets.  The 
^5-75%  duration  window  provides  good  agreement  between  observed  and  RVT 
predicted  peak  values. 

Estimation  of  esqiected  ground  motion  as  a  function  of  distance  and 
earldiquake  magnitude  is  fundamental  to  earthcpiake  hazard  assessment. 
The  general  problem  can  be  stated  as  follows:  given  an  earthquake  at  one 
location,  what  is  the  expected  ground  motion  at  other  locations?  Proper 
design  of  earthquake-resistant  structures  and  facdhties  requires  estimation  of 
the  ground  shaking,  typic^y  in  the  0.2  -  10  Hz  frequency  band.  In  addition, 
inelastic  dynamic  analysis  of  structures  requires  an  estimate  of  gignat  dura¬ 
tion. 

Ground  motion  at  a  particular  site  is  infiuenced  by  three  main  elements: 
source,  travel  path,  and  local  site  conditions.  Source  factors  include  size, 
depth,  stress  drop,  rupture  process  and  fault  geometry.  Travel  path  factors 
include  geometrical  attenuation,  dissipation  of  seismic  energy  due  to  anelas- 
tic  properties  of  the  earth,  and  focusing  and  scattering  of  elastic  waves  by  the 
three-dimensional  earth.  Local  site  factors  include  the  properties  of  the 
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uppermost  several  himdred  meters  of  rock  and  soil  and  the  effect  of  the  sur¬ 
face  topography  near  the  recording  site. 

The  study  of  ground  motion  attenuation  is  also  useful  for  the  inverse  pro¬ 
cess  of  using  groimd  motion  recordings  to  estimate  earthquake  source  param¬ 
eters,  such  as  magnitude  or  seismic  moment. 

Predictive  relationships  for  parameters  that  decrease  with  incresising 
distance  (such  as  peak  acceleration  and  peak  velocity)  are  often  referred  to  as 
attenuation  relationships.  A  number  of  attenuation  relationships  have  been 
proposed  for  western  and  eastern  North  America  in  the  past  two  decades 
(Joyner  and  Boore,  1981;  Boore  and  Joyner,  1991;  Boore,  1983;  Toro  and 
McGuire,  1987;  Atkinson  and  Boore,  1995;  Atkinson  and  Silva,  1997;  Camp¬ 
bell,  1981;  Campbell,  1985;  Campbell,  1997;  and  Sadigh  1997).  For  a  given 
region  these  relationships  may  differ  because  of  the  nature  of  the  data  sets, 
e.g.  three-component  accelerograms  or  regional  seismic  network  data. 

The  objective  of  this  study  is  to  characterize  ground  motion  observations 
in  Southern  California  from  the  TERRAscope  network  in  a  manner  consistent 
with  random  vibration  theory  (Boore,  1983),  which  requires  signal  duration 
and  amplitude  spectra  to  estimate  peak  motions.  Thus,  we  will  examine  the 
variation  of  Fourier  amplitude  spectra  with  distance  to  specify  a  geometrical 
spreading  function  and  a  frequency  dependent  Q  operator,  and  will  combine 
these  with  a  distance  and  frequency  dependent  duration  function  to  match 
peak  time  domain  observations. 

2.  Data  Set 

We  used  TERRAscope  seismograms  from  the  IRIS  Data  Management 
Center  to  characterize  distance  scaling  of  three-component,  0.25-5.0  Hz 
earthquake  ground  motions  in  the  distance  range  of  15-500  km  for  southern 
California.  The  data  set  consists  of  820  three-component  recordings  from  140 
earthquakes,  recorded  at  17  stations,  with  moment  magnitudes  between  3.1 
and  6.7.  The  maximum  frequency  of  5-Hz  was  constrained  by  the  20  sam¬ 
ples/sec  sampling  rate  of  the  broad-band  TERRAscope  data. 

Figure  1  shows  the  TERRAscope  stations  and  earthquakes  used  in  this 
study.  Table  1  gives  the  station  coordinates  and  a  description  of  the  site  geol¬ 
ogy  as  obtained  from  the  FDSN  station  book.  The  data  set  used  was  chosen 
to  have  good  overlapping  distance  sampling  by  stations  (Figure  2),  good  mag¬ 
nitude  coverage  (Figure  3),  and  independently  estimated  moment  magnitudes 
(Thio  and  Kanamoii,  1995).  The  TERRAscope  data  set  was  used  because  of 
its  broad  band,  the  broad  regional  coverage,  the  large  number  of  observa¬ 
tions  and  easy  accessibility.  The  broadband  channels  used  have  a  velocity 
sensitivity  that  is  flat  to  groimd  velocity  in  the  0.25  -  5.0  Hz  band. 
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Table  1 

_ TERRAscope  Station  Locatioiis  and  Geology 

Stations  Latitude  Longitude  Site  llype 

bar  32.680N  116.672W  ROCK,  Mesnznir  graniti7 

CALB  34.143N  118.627W 

DGR  33.650N  117.009W  ROCK,  Jurassic  metamorphic 

GLA  33.052N  114.827W 

GSC  35.303N  116.808W  ROCK,  Mesozoic  granitic 

ISA  35.663N  118.473W  ROCK,  Mesozoic  granitic 

MLAC  37.631N  118.834W  SOIL,  Quaternary  alluvium 

NEE  34.823N  114.596W  SOIL,  Quaternary  alluvium 

PAS  34.148N  118.172W  ROCK,  Cretaceous  quartz  diorite 

PFO  33.609N  116.455W  ROCK,  Cretaceous  granodiorite,  decomposed 

RPV  33.744N  118.404W  SOIL,  Quaternary  sedimentary  rock,  sand, rubble 

SBC  34.412N  119.713W  SOIL,  Quaternary  alluvium,  sand  Arid  gravel 

SMTC  32.949N  115.720W 

SNCC  33.248N  119.524W 

SVD  34.104N  117.097W  SOIL,  Qiiatemary  alluvium,  sand  and  gravel 

use  34.021N  118.287W  SOIL,  Quaternary  alluvium 

_ 34.567N  117.333W  SOIL,  Quaternary  alluvivun,  gravel,  sflT>d  and  gilf. 

Soil/rock  type  from  station  mformation  pages  at  IRIS  wHere  defined. 

Time  histories  were  obtained  from  the  IRIS  Data  Management  Center, 
corrected  for  instrument  response  to  yield  ground  velocity  in  m/s.  The  hori¬ 
zontal  components  were  rotated  into  radial  and  transverse  components.  The 
corrected  time  history  was  then  passed  through  a  program  that  filtered  the 
time  series,  determined  the  signal  duration,  Fourier  amplitude  spectrum,  and 
peak  motion.  Filtering  involved  passing  the  time  history  through  an  8-pole 
high-pass  followed  by  an  8-pole  low-pass  causal  Butterworth  filter  vdth  cor¬ 
ner  frequencies  of  0.707  and  1.414  ,  where  the  /*„ ’s  are  0.25,  0.33,  0.4, 

0.5,  1,  2,  3,  4,  and  5  Hz.  Both  Fourier  spectra  and  peak  time  domain  values 
are  used  in  the  regression  since  we  desire  a  self-consistent  model  linTrirtg  the 
three  observables. 

3.  Duration 

The  duration  of  ground  motion  can  have  a  strong  influence  on  earth¬ 
quake  damage.  A  motion  of  short  duration  may  not  produce  enough  load 
reversals  for  damaging  response  to  build  up  in  a  structure,  even  if  the  ampli¬ 
tude  of  the  motion  is  high.  On  the  other  hand,  a  motion  with  moderate  ampli¬ 
tude  but  long  duration  can  produce  enough  load  reversals  to  cause  substan¬ 
tial  damage  (Kramer,  1996). 

To  estimate  signal  duration,  integrated  square  filtered  velocity  is  used. 
Other  measinres  of  duration  have  been  defined  (Kramer,  1996),  but  we  wish  to 
use  one  consistent  with  random  vibration  theory.  Integration  starts  at  the  S- 
arrival  time  and  continues  into  the  coda.  The  duration  is  defined  as  the 
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window  within  which  the  integral  reached  5%  to  75%  of  its  ma-yimnm  This  is 
illustrated  in  Figure  4.  The  signal  duration  is  a  function  of  the  filter  center 
fi'equency  as  well  as  the  underl3ring  signal.  We  encoimtered  some  difficulties 
for  low  firequencies  because  of  microseism  noise,  which  introduced  a  linear 
trend  to  the  integral  at  large  time.  This  was  handled  by  examining  the  coda 
level  at  the  end  of  the  waveform  segment  to  limit  the  integration  window.  We 
also  noted  increasing  scatter  as  the  filter  firequencies  decreased  below  1.0  Hz. 

The  signal  within  these  5%  and  75%  limits  was  Fourier  transformed  and 
RMS  averaged  between  the  filter  corners  to  yield  the  Fourier  velocity  spectra 
in  m.  The  peak  velocity  value  of  the  filtered  time  history  following  the  S- 
arrival  time  was  also  saved.  The  data  set  thus  consists  of  observations  of 
duration,  Fourier  velocity  spectra  and  peak  values  of  filtered  velocities. 

A  consistent  definition  of  duration  is  required  because  an  objective  of  the 
study  is  to  have  a  set  of  observations  that  permits  the  use  of  random  vibra¬ 
tion  theory  as  described  by  Boore  (1983)  who  used  the  work  of  Cartwright 
and  Longuet-Higgins  (CLH)  (1956)  to  relate  Fourier  spectra  to  peak  motions 
through  signal  duration.  Equation  6.3  of  CLH  and  its  approximations  relate 
expected  peak  motion  to  RMS  estimates  by 

•^peak  ~  Vmax  -^RMSt  (1) 

where  Arms  is  estimated  using  ParsevaFs  theorem  and  the  filtered  Fourier 
spectrum  within  the  duration  window,  is  defined  by 

-1-00  ^ 

^max  ”  J  ^  d'H-  (2) 

— oo 

Here  N  is  the  number  of  maxima  to  be  exceeded,  and  q(ji)  is  the  cumulative 
probability  of  rj  exceeding  a  given  value  (CLH  5.1).  p(ji)  (CLH  1.19)  is  a  func¬ 
tion  of  the  spectral  moments;  and  the  number  of  peaks,  N,  depends  on  the 
moments  and  signal  duration.  Equation  6.2  of  CLH,  for  the  probability  den¬ 
sity  function  of  the  distribution  of  peaks, 

d 

pirimBx)  =  -j -  II  -  (3) 

is  used  to  examine  the  probability  of  not  exceeding  a  peak  value  by  providing 
confidence  bounds  on  Tjjaax' 

Figure  5  compares  the  ratio  of  observed  to  predicted  time  domain  peak 
values  at  each  filter  frequency.  For  each  observation  the  Fourier  velocity  spec¬ 
tra  for  the  time  window  is  used  with  the  duration  to  estimate  the  most  hkely 
peak  value  (Equation  2).  In  addition  the  5%  and  95%  bounds  of  the  cumula¬ 
tive  distribution  of  peaks  are  estimated  for  the  peak  using  Equation  3.  Fig¬ 
ure  5  aggregates  the  observations  according  to  these  bounds.  If  random 
vibration  theory  is  appropriate  and  if  the  duration  is  properly  defined,  then 
90%  of  the  observations  should  be  in  the  second  bin.  For  example,  at  2.00  Hz, 
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2033  of  the  2247  observations  are  in  the  second  bin.  Significant  deviations 
are  seen  at  the  low  and  high  filter  firequendes  for  data  with  low  values  of 
peak  motion.  Of  more  significance  is  the  observation  that  the  geometric  mpan 
of  the  ratio  of  observed  to  predicted  peeik  motions  is  0.97.  When  we  used  a 
longer  duration  based  on  the  5%  -  95%  integral  signal  squared,  this  ratio  dif¬ 
fered  significeintly  from  1.0.  Becaiise  of  this  exercise,  we  feel  that  we  have  an 
intemedly  consistent  data  set  of  a  duration  window,  Fourier  amplitude  spec¬ 
tra  and  peak  motions. 

Since  other  studies  (Atkinson  and  Boore,  1995)  support  a  distance-depen- 
dent  duration,  we  used  a  regression  model  consisting  of  station  terms  and 
piecewise  linear  segments  for  the  distance  dependence,  constrained  to  be  zero 
at  zero  distance  and  to  be  smooth.  Station  terms  were  necessary  to  account 
for  the  variability  of  durations  at  low  filter  frequencies.  We  found  significant 
fi'equency  dependence  in  the  quahty  of  the  duration  data  set  and  trends.  Fig¬ 
ures  6  and  7  show  the  raw  three-component  duration  data  for  the  0.25  Hz 
and  1.00  Hz  filter  frequencies  as  well  as  the  residuals  based  on  the  model. 
The  low  frequency  duration  data  exhibit  much  scatter  up  to  0.5  Hz.  In  addi¬ 
tion,  the  station  residuals  show  interesting  patterns. 

Table  2  presents  the  station  terms  at  frequencies  of  0.25  Hz  and  1.0  Hz. 
To  perform  the  regression,  the  station  correction  at  Goldstone  (GSC)  was  con¬ 
strained  to  be  zero.  The  Sigma  value  is  the  standard  error  of  the  mean  Sta¬ 
tion  Term  and  not  an  indication  of  distribution  of  residuals  themselves.  For 
stations  with  more  than  90  observations,  we  note  for  example,  that  Santa 
Barbara  (SBC)  has  a  longer  duration  than  Pasadena  (PAS)  at  at  0.25  Hz.  The 
vanabihty  in  station  terms  may  be  an  artifact  of  the  effect  of  microseism 
noise  on  the  automatic  determination,  but  may  also  be  real,  refiecting  rever¬ 
beration  of  these  low  frequency  waves  caused  by  local,  perhaps  shallow,  3-D 
structure.  The  standard  error  of  fit  of  53.5  sec  is  high  for  the  0.25  Hz  data, 
perhaps  because  no  data  screening  was  performed.  Figure  6  seems  to  show  a 
tendency  for  the  scatter  in  residuals  to  increase  with  distance.  The  regres¬ 
sion  model  of  station  terms  plus  a  distance  effect  may  be  inappropriate  in  the 
sense  that  the  station  terms  are  assumed  independent  of  distance. 

The  conclusion  is  that  station  terms  are  required  and  that  both  the 
degree  of  distance  dependence  and  scatter  decrease  with  increasing  fre¬ 
quency.  In  a  comparison  of  durations  at  all  filter  frequencies.  Figure  8  shows 
that  the  low  frequency  signal  components  have  longer  durations  than  high 
frequency  components.  The  durations  at  different  distances  for  different  fre¬ 
quencies  are  Ksted  on  Table  3. 
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Table  2 

Station  Residuals  fix)m  Duration  Regression  Analysis 

Station _ Station  Tterm _ Sigma  NOBS 

0.25  Hz  Variance  Reduction  17.5%  cr  53.5  sec 


BAR 

25.3 

10.8 

31 

CALB 

24.7 

11.9 

21 

DGR 

70.4 

15.5 

13 

GLA 

29.5 

20.9 

7 

GSC 

0.0 

0.0 

119 

ISA 

18.1 

6.4 

129 

MLAC 

8.6 

12.6 

21 

PAS 

25.3 

5.2 

178 

PFO 

24.5 

5.5 

133 

RPV 

61.3 

8.8 

43 

SBC 

59.7 

6.9 

94 

SNCC 

0.1 

31.5 

3 

SMTC 

63.9 

19.1 

9 

SVD 

18.3 

7.8 

58 

use 

57.3 

9.1 

39 

VTV 

12.3 

8.5 

48 

1.00  Hz  Variance  Reduction  53.5%  a  6.7  sec 


BAR 

3.9 

1.0 

59 

CALB 

9.8 

1.0 

46 

DGR 

2.1 

1.8 

15 

GLA 

-3.6 

2.0 

12 

GSC 

0.0 

0.0 

234 

ISA 

-1.1 

0.6 

281 

MLAC 

7.4 

1.2 

37 

NEE 

21.8 

1.4 

24 

PAS 

0.0 

0.5 

333 

PFO 

3.8 

0.6 

255 

RPV 

8.3 

0.8 

79 

SBC 

11.7 

0.7 

174 

SNCC 

16.1 

2.7 

6 

SMTC 

41.5 

2.1 

12 

SVD 

2.4 

0.7 

113 

use 

17.1 

0.8 

77 

VTV 

0.6 

0.8 

86 

4.  Ground  Motion  Regression 


The  regression  analysis  for  peak  motion  and  Fourier  velocity  spectra 
uses  the  same  simple  model.  Assuming  a  multiplicative  effect  of  source, 
propagation,  and  site,  the  logarithm  of  the  ground  motion  parameter  is  mod¬ 
eled  as  the  separable  effects  of  source,  site,  and  propagation: 
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Tables 


Distance  _ Durations  for  different  Frequencies  (Hz) 


(km) 

0.25 

0.33 

0.4 

0.5 

1.0 

2.0 

3.0 

4.0 

5.0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

40 

12.47 

9.81 

12.30 

13.12 

9.16 

5.90 

4.78 

4.09 

3.72 

80 

24.63 

19.37 

22.80 

20.23 

11.55 

7.86 

6.60 

5.67 

5.04 

120 

35.75 

28.21 

28.34 

21.01 

10.61 

8.06 

6.64 

5.71 

5.09 

160 

45.21 

35.99 

28.96 

21.26 

12.33 

9.08 

7.41 

6.36 

5.61 

200 

53.20 

42.90 

36.39 

28.49- 

14.82 

10.96 

8.99 

7.91 

7.11 

250 

60.10 

49.05 

37.71 

28.62 

17.10 

13.10 

11.09 

9.66 

8.90 

300 

66.74 

54.92 

33.74 

24.65 

19.20 

14.96 

13.24 

11.31 

10.72 

PEAK = log  A = SRC(.f) + SlTEif)  +  D(r,  f)  (4) 

Here  D(r)  is  the  propagation  term,  expressing  the  combined  effect  of  geo¬ 
metrical  spreading  and  attenuation.  While  simple  in  appearance,  true  sepa- 
rabihty  is  impossible  because  of  tradeofife.  Care  must  be  taken  to  select  sites 
which  observe  earthquakes  at  many  distances  so  that  observations  from  sites 
overlap  to  avoid  an  undesirable  trade-off  between  an  event  source  term  and 
the  distance  function.  Constraints  must  also  be  applied,  such  as  forcing 
D(r)  =  0  at  some  reference  distance,  such  as  Dir  ref)  -  0.  The  reference  dis¬ 
tance  Tref  used  is  such  that  errors  in  source  depth  make  little  difference  on 
hypocentral  distance  and  that  regional  variations  in  high  amplitude  Moho 
reflected  signals  are  avoided.  Additionally,  some  or  all  of  the  site  terms  must 

n 

be  constrained,  such  as  X  SITEi = 0 ,  which  has  the  side  effect  of  forcing  com- 

M) 

mon  site  effects  into  all  source  terms.  This  permits  the  determination  of  a 
relative  site  response  only.  Any  common  site  term  is  mapped  into  the  source 
term.  For  this  study,  this  constraint  was  placed  only  on  the  horizontal  compo¬ 
nents  -  the  vertical  components  were  free  to  move  relative  to  the  horizontal. 
Since  the  initial  computations  show  that  the  SITE  terms  for  the  horizontal 
components  at  GSC  were  small,  we  further  constrained  each  of  the  two  hori¬ 
zontal  component  SITE  terms  for  GSC  to  be  zero. 

To  emphasize  the  fact  that  we  are  not  defining  the  earthquake  source  on 
the  basis  of  our  data  set,  but  rather  parameterizing  the  observed  motions, 
equation  (4)  is  rewritten  as 

PEAK  =  log  A = Eif,  Vref ) + S(/)  +  D(r,  /) ,  (5) 

where  E  represents  the  excitation  of  the  ground  motion  at  the  reference  dis¬ 
tance  and  S  is  the  site  term.  Later  we  will  compare  predicted  levels  of 
ground  motion  at  the  reference  distance. 

Investigations  by  Atkinson  and  Boore  (1995)  in  eastern  North  America 
and  by  Atkinson  and  Silva  (1997)  in  Cahfomia,  indicate  that  the  distance 
dependence  does  not  have  a  simple  functional  form.  However,  the  geometri¬ 
cal  spreading  may  be  expressed  as  a  piecewise  continuous  linear  relation. 
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Anderson  and  Lei  (1994),  Savage  (1995)  and  Harmsen  (1997)  extended  this 
concept  to  permit  many  linear  segments  to  be  tied  together  with  a  smooth¬ 
ness  constraint.  This  means  that  definition  of  geometrical  spreading  and  fi*e- 
quency  dependent  Q  is  deferred  to  a  later  stage.  At  a  given  firequency  we 
model  D(r)  as  a  piecewise  linear  function  defined  by 

D(r)  =  tLi(r)Di,  (6) 

1=0 

and  where  Lj  (r)  is  a  linear  interpolation  function,  and  the  D,-  are  node  values 
such  that  D(ri)  =  Dj.  A  smoothness  constraint  can  be  applied  by  requiring 

Di_i  -  2Di  +  Di^i  =  0.  (7) 

Minimum  roughness  is  easily  incorporated  if  the  nodes  at  distances  are 
evenly  spaced. 

In  this  study,  a  coda  normalization  (AM,  1980  ;  Frankel  et  al,  1990)  was 
used  to  provide  an  initial  estimate  of  the  distance  term  Z)(r),  by  removing  the 
effects  of  somce  and  site  jfrom  the  logarithm  of  peak  groimd  motion,  A.  The 
initial  D(r)  was  then  used  in  an  interactive,  damped  least-squares  regression 
to  estimate  excitation,  site  and  distance  functions.  Alternatively,  a  high  qual¬ 
ity  data  set  such  as  this  one  could  start  with  jD(r)=logio(r„//r)  as  an  initial 
estimate.  The  distance  function  is  parameterized  as  a  piecewise  linear  func¬ 
tion  with  17  nodes  between  10  and  500  km. 

For  our  regression,  we  chose  r^ef  =  40  km.  For  each  filter  firequency  there 
were  140  event  terms,  51  site  terms  for  the  17  three-component  stations,  and 
13  nodes  in  the  distance  function.  Table  4  gives  the  number  of  observations, 
standard  deviation  of  residuals  and  standard  error  of  the  mean  residual  for 
the  Fourier  and  time  domain  regressions.  The  mean  residual  is  zero. 

’Ibble4 

Regression  Error  Analysis 
Fourier  Domain  Time  Domain 

fn  (Hz)  Nqbs  ^OBS  ^ 


0.25 

2204 

0.299 

0.006 

2204 

0.153 

0.003 

0.33 

2284 

0.162 

0.003 

2284 

0.155 

0.003 

0.40 

2349 

0.115 

0.002 

2349 

0.158 

0.003 

0.50 

2375 

0.117 

0.002 

2375 

0.158 

0.003 

1.00 

2408 

0.125 

0.003 

2408 

0.164 

0.003 

2.00 

2403 

0.135 

0.003 

2403 

0.180 

0.004 

3.00 

2384 

0.142 

0.003 

2384 

0.185 

0.004 

4.00 

2354 

0.150 

0.003 

2354 

0.187 

0.004 

5.00 

2317 

0.161 

0.003 

2317 

0.193 

0.004 

Figures  9  and  10  show  the  D(r)  functions  for  the  Fourier  and  time 
domain  data  sets,  respectively,  as  a  function  of  frequency.  To  enhance  presen¬ 
tation,  these  figures  show  the  deviation  firom  a  r~^  trend.  For  clarity,  error 
bars  are  not  plotted.  Subtle  differences  are  seen  at  larger  distances,  which 
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are  interpreted  as  the  effect  of  duration  on  the  time  domain  observations. 
Figures  11  and  12  show  the  distribution  of  regression  residuals  as  a  function 
of  distance  for  the  two  data  sets  at  selected  frequencies.  This  plot  is  used  to 
check  the  appropriateness  of  the  distance  nodes  used  to  determine  Z)(r)  since 
distance  trends  would  be  apparent. 

The  site  terms  are  similar  for  the  time  and  frequenty  domain  data  sets. 
Figures  13  and  14  give  the  site  terms  for  the  Fourier  velocity  spectra  data 
set.  In  general,  the  site  terms  for  the  radial  and  transverse  components  over¬ 
lay.  The  horizontal  components  site  terms  exceed  the  vertical  terms  and  the 
frequency  dependence  varies  from  station  to  station.  Because  of  the  con¬ 
straint  apphed  to  the  horizontal  components  site  term  at  GSC,  we  may  inter¬ 
pret  the  excitation  terms  of  equation  (4)  as  being  the  horizontal  ground 
motion  level  for  an  earthquake  40  km  from  GSC.  This  distinction  is  impor¬ 
tant  if  one  wishes  to  relate  the  excitation  and  site  terms  to  strong  motion 
data  sets,  which  use  data  from  a  different  network. 

There  is  much  variation  in  the  site  terms  among  the  stations  shown  in 
Figures  13  and  14.  The  extremes  are  bounded  by  [-0.7,  0.6].  The  horizontal  to 
vertical  ratio  is  sKghtly  less  variable.  This  logarithmic  ratio  of  horizontal  to 
vertical  component  site  terms,  and  mean  logarithmic  ratios  for  soil,  rock,  and 
all  site  types  are  shown  on  Figure  15  and  are  listed  in  Table  5,  for  the  Fourier 
velocity  spectra.  The  site  type  is  that  given  in  Table  1.  There  is  little  differ¬ 
ence  between  the  estimates  for  the  two  domains.  The  H/Z  ratio  mean  varies 
for  the  different  site  types  with  a  mean  value  corresponding  to  about  V2. 

Table  5 


Mean  Logarithmic  Ratio  of  Horizontal  lb  Vertical  Site  Terms 


/.(Hz) 

Rock 

All 

Soil 

0.25 

0.138 

0.236 

0.33 

0.131 

0.237 

0.40 

0.127 

0.50 

0.125 

1.00 

0.171 

0.225 

2.00 

0.258 

3.00 

0.237 

4.00 

0.169 

0.222 

5.00 

0.175 

0.179 

5,  Ground  Motion  Parameterization 

The  next  step  in  the  analysis  is  to  define  geometrical  spreading  and 
anelastic  attenuation  functions  to  describe  the  Fourier  velocity  spectra  Dir) 
term.  TMs  is  done  by  assuming  Q(f)=Q()f''  and  a  simple  piecewise  linear 
geometrical  spreading  function.  The  frmctional  form  g(r)  for  geometrical 
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spreading,  is 


(r/r 

\*  f  •  cross/ 


r<r, 
r  >  r 


cross 


cross 


The  effect  of  anelastic  attenuation  is  to  reduce  amplitude  with  distance  by  a 
factor  of  ex^-7cfrlQif)p),  where  >?  =  3. 5  km/sec. 

For  a  set  of  the  parameters  Qq,  rj  and  r gross,  theoretical  Fourier  spectra 
were  estimated  at  each  of  the  distances  used  in  the  regression  for  Dir\  the 
results  were  normalized  to  the  reference  distance  of  40  km,  and  a  logarithm 
taken  for  direct  comparison  with  the  regression  results.  A  search  is  made 
through  the  parameter  space  to  find  the  rmniTnnm  of 


i 

X(l/<Ti)2 

i 

The  <7,  are  those  resydting  from  the  regression.  As  a  check  we  also  compared 
the  Dir)  values  obtained  irom  the  time  domain  regression  by  making  random 
vibration  theory  estimates  of  peak  filtered  groimd  velocities  using  the 
observed  durations  and  random  vibration  theory  (Boore,  1983).  We  accounted 
for  the  time  domain  response  of  the  filters  in  the  manner  of  Boore  and  Joyner 
(1984)  for  lightly  damped  single  degree  of  fireedom  osdllators  by  stating  that 
the  R]\^  dxiration  is  the  smn  of  the  soiurce  duration,  propagation  duration, 
and  twice  the  filter  period.  In  addition  the  duration  used  for  determining  the 
number  of  random  peaks  is  the  sum  of  the  source  duration,  the  propagation 
duration  and  the  filter  period. 

This  exercise  yielded  =40  km,  Qo  =  180  and  7  =  0.45.  Figures  16 
and  17  present  the  residuals  for  the  Fourier  velocity  and  timp  domain  Dir) 
fimctions,  respectively.  The  comparison  shows  the  values  for  the  distance 
range  of  20  -  500  km,  because  of  the  paucity  of  data  in  the  10  -  20  km  range. 
The  fit  of  the  time  domain  data  is  an  indication  of  the  appropriateness  of  the 
duration  function  used,  when  combined  with  the  firequency  domain  model. 

Previous  work  by  Atkinson  and  Silva  (1997)  resulted  in  Qo  =  204, 
»7=0.56,  and  a  trilinear  geometric  spreading  fimction  with  exponents  of -1  0 
and  172  for  r  ^  50  km,  50  ^  r  <  170  km,  and  r  >  170  km,  respectively.  A  com- 
panson  of  the  difference  between  our  observed  Fourier  velocity  Dir)  and  that 
predicted  by  their  values  is  given  in  Figure  18.  Their  data  set  has  relatively 
few  observations  beyond  100  km  and  is  dominated,  by  the  June  28,  1992, 
Mjf,  -  7. 3,  at  large  distances.  We  have  better  agreement  with  the  regressions 
of  Harmsen  (1987)  who  used  TERRAscope  data  in  the  10  -  150  km  range. 
The  relation  for  Q(f)  obtained  in  oiu:  study  is  more  in  agreement  with  the 
regional  attenuation  relationship  for  southern  California  by  Benz  et  al 
(1997),  which  gives  Q  if)  =  187!?  However,  their  data  set  includes 

both  southern  and  central  California  for  a  distance  range  of  150 <r< 700  Irm 
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As  a  final  comparison,  we  compared  the  observed  exdtation  levels  of  the 
horizontal  motion  for  both  data  sets  at  40  km  to  model-based  predictions  as  a 
function  of  moment  magnitude.  Since  there  are  few  data  for  Mw  >  5. 5,  we 
did  not  feel  confident  in  defining  a  somce  exdtation  model,  which  would  also 
have  suffered  by  the  extrapolation  of  our  observations  to  a  nominal  1  km  dis¬ 
tance.  Instead  we  compared  our  exdtation  levels  at  40  km  to  that  derived 
fi:om  two  models  in  literature,  one  by  Atkinson  and  Boore  (1998)  and  the 
other  by  Boore  and  Joyner  (1997).  The  model  for  the  horizontal  site  accelera¬ 
tion  spectrum  is 

A(f)  =  C(.2^ffMoS(f)g(r)expi-^rlWfW(.f) 

where 

C  =  (0. 55)(0. 707)(2. 0)/47cpf , 

and  V(f)  is  a  site  amphfication  term.  The  fimctional  form  for  Sif)  is 

1  —  f  € 

and  the  /«,  and  e  are  given  by  Atkinson  and  Boore  (1998)  2ind  repeated  in 
Table  6. 

Table  6 

Source  Parameters 


Atkinson  and  Boore  (1998) 

p  2.8 

p  3.5 

fa  log/“a=2.181-0.496M  M>4.8 
=2.617-0.500M  M<4.8 
fb  logA  =  1.308-0.277M  M>4.8 
=2.617-0.500M  M<4.8 
e  logs  =  3.223-0.670M  M^4.8 
=0.0  M<4.8 

Boore  and  Joyner  (1997) 

p  2.8 

p  3.5 
Act  70  bars 

/•«=4.9x10V(A<t/Mo)'^" 

fh  fb-fa 

e  g=1.0 _ . 

We  use  propagation  durations  of  12.5  sec  at  0.25  Hz,  9.2  sec  at  1.0  Hz,  and 
3.7  sec  at  5  Hz  at  40  km.  The  source  duration  is  0. 5//‘a.  We  use  a  x*  of  0.040 
sec  and  0.035  sec  for  the  Atkinson  and  Boore  (1998)  and  Boore  and  Joyner 
(1997)  models,  respectively.  The  Vif)  term  is  taken  fi:om  Atkinson  and  Silva 
(1997)  and  is  based  on  that  in  Boore  (1986). 


12 


Figures  19  and  20  compare  the  excitations  at  40  km  from  our  regressions 
with  the  predictions  based  on  the  two  somce  models  combined  with  our  geo¬ 
metrical  spreading,  anelastic  attenuation  and  duration  functions.  To  empha¬ 
size  differences,  E-vS-M^  is  plotted.  Om:  results  share  the  trends  of  the 
two  models  as  a  function  of  moment  magnitude.  Both  models  imderpredict 
low  frequency  levels  at  small  moments.  The  Boore  and  Joyner  constant  stress 
drop  (1997)  model  seems  to  fit  the  high  frequency  data  better  than  the  Atkin¬ 
son  and  Boore  (1998)  model,  but  recall  that  these  figures  are  for  expected 
motions  at  GSC. 

Finally,  the  results  of  this  study  indicate  that  the  modem  TERRAscope 
data  set  agrees  with  strong  motion  measmrements  in  estimating  distance 
dependence,  which  has  been  also  noted  by  Harmsen  (1997).  A  direct  compari¬ 
son  with  strong  motion  data  sets  requires  a  joint  inversion  of  both  data  sets. 

6.  Conclusions 

Using  ground  velocities  from  the  three-component  recordings  of  the  TER¬ 
RAscope,  we  have  characterized  the  groimd  velocity  distance  scaling  in  the 
range  of  15  -  500  km  for  Southern  California.  Our  groimd  motion  scaling  is 
similar  to  those  of  other  studies  whidi  have  used  TERRAscope  data  (Benz  et 
al.,  1997;  Harmsen,  1997).  The  major  differences  of  our  study  are  the  use  of 
Fourier  velocity  spectra,  signal  duration,  and  peak  value  of  filtered  time 
domain  velocities. 

The  results  of  this  study  indicate  that  low  frequency  signal  components 
have  longer  durations  than  high  frequency  components  and  increase  signifi¬ 
cantly  with  distance.  The  low  frequency  duration  data  exhibit  much  scatter 
up  to  0.5  Hz.  The  duration  results  indicate  that  station  terms  are  required 
and  that  both  the  degree  of  distance  dependence  and  scatter  increase  with 
decreasing  frequency.  The  variability  in  the  station  terms  may  be  due  to  the 
effect  of  micxoseism  noise  or  due  to  reverberation  of  low  frequency  waves  by 
local  3-D  structure.  These  long  durations  may  be  important  for  assessing  the 
response  of  stmctures  with  low  natural  frequencies. 

The  distance  fimctions  obtained  from  the  Fourier  and  time  domain  data 
sets  are  similar  in  trend.  Subtle  differences  at  larger  distances  can  be  inter¬ 
preted  as  the  effect  of  duration  on  the  time  domain  observations.  The  site 
terms  are  similar  for  the  time  and  frequency  domain  data  sets.  In  general, 
the  site  terms  for  tiie  radial  and  transverse  components  are  quite  similar. 
The  horizontal  motions  exceed  the  vertical  and  the  frequency  dependence  of 
this  ratio  varies  from  station  to  station.  The  ratio  of  horizontal  to  vertical 
component  site  terms  is  about  V2  for  all  frequencies.  However,  this  ratio  is 
close  to  one  for  rock  sites  but  higher  for  soil  sites  alone.  There  is  much  vari¬ 
ability  in  the  site  terms  among  stations. 
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Random  Vibration  Theory  (RVT)  has  been  used  to  model  the  observed 
peak  ground  motion.  Peak  velocities  are  controlled  by  a  combination  of  dura¬ 
tion,  geometrical  spreading,  and  anelastic  attenuation.  The  results  of  RVT 
appHcation  indicate  that  our  data  fit  well  by  bilinear  geometrical  spreading 
e^onents  of  -1  for  r  <  40  km  and  -1/2  for  r  >  40  km.  The  frequency  depen¬ 
dent  quahfy  factor  is  Q  (/)  =  180  for  the  combined  three-component  data 
sets.  Omr  observations  show  the  effect  of  duration  as  a  factor  relating  the 
Fourier  amplitude  spectra  to  the  peak  amplitude  seeing. 

The  comparison  of  the  excitation  at  40  km  from  our  regressions  ivith  the 
predictions  based  on  the  two  models,  one  by  Atkinson  and  Boore  (1998)  and 
the  other  by  Boore  and  Jo3mer  (1997),  indicate  that  our  results  share  the 
trends  of  the  two  models  as  a  function  of  moment  magnitude.  Both  models 
overpredict  high  frequency  levels.  This  is  due  to  our  use  of  the  V(f)  function 
developed  for  strong  motion  sites  and  the  fact  that  the  TERRAscope  data  set 
has  not  been  used  together  with  strong  motion  data  in  an  inversion.  Because 
of  the  smaller  number  of  large  events,  we  cannot  distinguish  between  the 
Atkinson  and  Boore  (1998)  and  the  Boore  and  Joyner  (1997)  models,  even 
though  the  Boore  and  Joyner  (1997)  model  seems  slightly  better. 

A  comparison  of  the  Fourier  domain  excitation  terms  (Fig.  19)  to  the  tiVna 
domain  terms  (Fig.  20)  shows  the  same  relative  pattern  between  the  data  anH 
predictions.  This  supports  the  internal  consistency  of  our  pflrflTnptAri7:f^t■■ior^  of 
the  data  in  terms  of  anelastic  attenuation,  geometric^  spreading,  and  dura¬ 
tion,  since  the  time  domain  values  use  a  duration  that  incx’eases  with 
moment  magnitude. 

Finally,  the  results  of  this  study  indicate  that  the  modem  TERRAscope 
data  complement  strong  motion  measurements  in  estimating  distance 
dependence  (Figure  18  at  distances  less  than  80  km). 
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Fig.  1.  Locations  of  TERRAscope  stations  used  (solid  triangles)  and  earth¬ 
quakes  (open  circles).  The  size  of  the  circle  is  an  indication  of  relative 
earthquake  magnitudes. 
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Fig.  2.  Distance  distribution  of  data  set  arranged  by  station. 


Fig.  3.  Distance  distribution  of  data  set  arranged  by  moment  magnitude. 
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Fig.  4.  Illustration  of  method  of  estimating  duration.  The  lower  trace  is  the 
velocity  time  history  filtered  at  1.0  Hz.  The  upper  trace  is  the  integrated 
square  velocity.  The  dirration  of  14  s  is  the  time  interval  between  the  0.05 
and  0.75  ordinate  (small  squares  and  heavy  vertical  bars  in  upper  figure). 
The  origin,  P-  and  S-arrival  times  (IPUO  and  ISUO,  respectively)  from  the 
unfiltered  time  history  are  indicated.  The  group  delay  of  the  filter  is 
apparent  by  the  shift  of  the  P-  and  S-arrivals  from  the  picked  time. 
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PREDICTED  PEAK  PREDICTED  PEAK 


Fig.  5.  Comparison  of  ratio  of  observed  to  predicted  peak  amplitudes  as  a 
function  of  predicted  peaks  for  the  entire  data  set  at  the  different  filter 
firequendes.  Each  figure  indicates  the  filter  fi:equen(y,  the  total  number  of 
observations  separated  as  to  whether  tiie  observed  peak  was  in  the  0-5%, 
5-95%  or  in  the  95-100%  bounds  of  the  predicted  peak,  the  distance  range’ 
and  the  geometric  mean  ratio.  ’ 
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Fig.  6.  Duration  data  and  regression  line  (top)  and  regression  residuals  as  a 
function  of  distance  for  a  filter  frequency  of  0.25  Hz. 
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Hypocentral  Distance  (km) 


Fig.  7.  Dxxration  data  and  regression  Kne  (top)  and  regression  residuals 
function  of  distance  for  a  filter  firequency  of  1.00  Hz. 
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Fig.  8.  Distance  dependence  of  duration  for  0.25,  0.40,  1.0,  3.0  and  5.0  Hz  fil¬ 
tered  data. 
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HYPOCENTRAL  DISTANCE  (km) 

Pig.  9.  Frequency  dependence  of  Fourier  velocity  spectra  D(r)  corrected  for 
spreading. 
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Fig.  10.  Frequency  dependence  of  filtered  time  domain  D{r)  corrected  for  BT^ 
spreading. 
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Fig.  13.  Site  terms  for  Fourier  velocity  spectra  as  a  function  of  frequency.  The 
Z  component  is  indicated  by  the  dark,  solid  line;  the  R  component  is  indi¬ 
cated  by  the  medium  shade,  short-dashed  line  ;  the  T  component  is  indi¬ 
cated  by  the  light,  long-dashed  line.  Regression  error  bars  are  indicated. 
Adjacent  to  the  station  name  is  the  number  of  vertical  component  observa¬ 
tions  at  l.OHz;  this  permits  a  visual  weighting  of  the  importance  of  any 
station  site  term. 
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Pig.  14.  Site  terms  for  Fourier  velocity  spectra  (cont’d). 
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Fig.  15.  Observed  and  mean  HIZ  ratio  from  the  station  terms  of  Fig.s  13  and 
14. 
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HYPOCENTRAL  DISTANCE  (km) 

Fig.  16.  Residuals  of  the  model  fit  to  the  Fourier  velocity  distance  function. 


30 


o 

w 


HYPOCENTRAL  DISTANCE  (km) 

Fig.  17.  Residuals  of  the  model  fit  to  the  filtered  time-domain  velocity  dis¬ 
tance  function. 
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HYPOCENTRAL  DISTANCE  (km) 


PiS-  18,  Residuals  of  tlie  model  fit  to  the  Fourier  velocity  distance  function 
using  the  Atkinson  and  Silva  (1997)  Qif)  and  geometrical  spreading  func¬ 
tion. 
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Pig.  19.  Excitation  of  Fourier  velocity  spectra  at  40  km.  The  solid  curve  is  the 
prediction  of  the  Atkinson  and  Boore  (1998)  source  model;  the  deished 
curve  is  that  of  the  Boore  and  Joyner  (1997)  source  model.  E  is  the  log^o  of 
the  Fourier  velocity  spectra  in  m. 
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Fig.  20.  Excitation  of  peak  filtered  time  domain  velocities  at  40  km.  The  solid 
curve  is  the  prediction  of  the  Atkinson  and  Boore  (1998)  source  model;  the 
dashed  curve  is  that  of  the  Boore  and  Jo3mer  (1997)  source  model.  E  is  the 
logio  of  the  filtered  velocity  spectra  in  ml  sec. 
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